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ТРЕНИЕ И ИЗНОС 


УДК 681.51. 
В.Л.ЗАКОВОРОТНЫЙ, А.Д.ЛУКЬЯНОВ, В.Е. СТУПИН 


УСЛОВИЯ ПОТЕРИ УСТОЙЧИВОСТИ ДВИЖЕНИЯ ы 
МЕХАНИЧЕСКОЙ СИСТЕМЫ, ВЗАИМОДЕИСТВУЮЩЕИ 
С ТРИБОСРЕДОИ 


Рассматриваются условия потери устойчивости движения механической системы, 
взаимодействующей со средой, формируемой в трибосопряжении. Показаны усло- 
вия возникновения и роль диссипативных, гироскопических, потенциальных и цир- 
куляционных сил, формируемых в узле трения естественным образом. 

Ключевые слова: трение, устойчивость движения, управление. 


Введение. При проектировании машин и механизмов часто необходимо 
решать вопросы устойчивости движения с учетом существующих трибоузлов. 
В зоне, находящейся между контактируемыми поверхностями трибоузла, фор- 
мируется диссипативная среда, образующая динамическую связь [1-3]. Дина- 
мические свойства такой системы раскрываются на основе «базовой» дина- 
мической модели трибосистемы [4], рассматривающей смещение упруго- 
диссипативно подвешенного индентора по отношению к абсолютно жёстко- 
му основанию (контробразцу) большой массы (рис.1). Индентор представ- 
ляется в виде сосредоточенной массы. Он движется с медленно меняющей- 
ся скоростью И (Ё), задаваемой управлением от внешнего источника. Урав- 
нение движения такой системы можно представить в виде 


тХ + ВХ + с = Е(Х, Х.Г] + Е' (0+0, (1) 


ге Х={Х,Х,\' - вектор координат системы. Координаты вектора 


состояния показывают текущие положения поверхности индентора в 
трибосреде (формируемой в процессе относительного скольжения между 

 ЕЖ.Х = Ех, ,Е ХХ 
контактируемыми поверхностями); - 
вектор-функция, описывающая зависимость сил контактного 
взаимодействия от координат состояния. Её свойства определяются 
формируемой в процессе трения переходной областью между индентором 


и контробразцом, названной трибосредой; ЕР” (1) - силы, не объяснимые в 


координатах состояния системы, которые можно интерпретировать как 
неуправляемый СИЛОВОЙ шум, возмущающий стационарные движения 


системы; т,й,с - матрицы размерности 2 2 соответственно обобщенных 

















масс, коэффициентов демпфирования и жесткости; | “< | - вектор 
внешних силовых воздействий. Таким образом, внешние условия 
функционирования узла трения задаются | < |и И(@) - медленными 











функциями времени. 
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Е=4Е (ХЕ (ХУ 


Рис.1. Схема «базовой» динамической модели узла трения 


Ограничимся наиболее важным для технических приложений 
случаем, когда стационарная траектория в пространстве состояния есть 

















точча Х’(Й=Х’, что — соответствует с = = соп51 и 
(0 ЕТ = сопт5{. Возможность такого рассмотрения реальной системы 
определяется тем, что |_“® | ИС) - медленные функции времени, и 











изменения координат Х”(1] при условии асимптотической устойчивости 
системы также будут описываться медленной функцией времени. Тогда, 


























учитывая, что ао =—- “емо |. где = - Вариации относительно 
стационарной точки | х | получаем: 
тХ+ ВХ сх ХИ, хх] + Е" (0, (2) 


где © Йа = жж [ХИ] - нелинейные функции 


в вариациях относительно точки равновесия, причем ф (Х”,И,0,0) = 0. 
Т Т 


Кроме этого справедливо Х,4Ё-> 0 и х,4Ё-> О при достаточно 


0 0 
большом 7’, в противном случае Х’- . Для стационарных 
периодических траекторий с периодом 7’ также выполняются условия: 
Т Г. 
х@=0и Хх, -= 0, 
о 0 


в противном случае после каждого периода колебаний формируется 
дополнительное смещение точки равновесия. Случай реверсивного трения 
не рассматривается. 

Для автономной системы: 


тх+ Их с х= 0, (3) 


где № = 1,,-099,/9х, , <= с,,-960,/9х, ‚ 5,К=12 устойчивость 











траектории | х`х | определяется свойствами матриц Я; и с... 
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Подчеркнем, что 90, 9х, и д0,/0х, зависят от траекторий | > |и 





динамической характеристики @ (№,,х,х) в вариациях относительно 


этой траектории. Поэтому изменение Х’”, обусловленное изменениями 
внешних условий С/ или И, вызывает вариации Я; и с, и 
проявляется в смещении корней характеристического полинома. Все может 
привести к принципиальному изменению свойств системы (3) вплоть до 
потери точкой равновесия Х’* устойчивости. 

Исходная динамическая система без учета трибосреды задается 
матрицами Ее} Я [.] (левые части в (2)), которые имеют постоянные 
параметры и являются положительно определенными и симметричными. 
Без трения исходная система имеет асимптотически устойчивую точку 
равновесия. Диссипативную функцию введем в форме Релея. 

При переходе к (3) матрицы и ‹ | уже не обладают 
свойством симметрии, так как в общем случае 90,/9х, 90,/9х, и 
90,/9х, 960,/9х,, 5К=1,2. Поэтому свойства (3) принципиально 
меняются. Матрицы | и Ш можно представить как 















































| ее не” р р; = И + И , (4) 
ы Т ы т 
где <= 1/2 с.+(с) , №?=12 и +(№} - симметричные части матриц 
жёсткости и диссипации, отвечающие за потенциальные свойства 
(к) Т 
= 2 в -[с (к) _ е т 
системы; ^ Рае, ; В” = УЗИ - №} |. кососимметричные 
матрицы, имеющие структуру 
: (К) . (К) 
(к) _ 0; = Си (к) _ 0; д й; 1,5 
а 0 и № = р. 0 
С; 51552 У ,5 2 
Характерной особенностью сил, формируемых матрицей | * | 
является то, что их работа по замкнутому контуру, содержащем 
стационарную точку, равна нулю. Силы, формируемые матрицей | № | 











если она является положительно определенной, всегда направлены против 
движения и совершают не нулевую работу. При этом силы, зависящие от 






















































































матрицы | =” | не совершая работу, могут способствовать стабилизации 
или раскачиванию точки равновесия. 

Что касается матриц |“ | и < | то силы, формируемые 
элементами матрицы | <’ совершают работу при движении по 
замкнутому контуру, а формируемые элементами матрицы | <’ | не 
совершают работу по замкнутому контуру. Более того, если матрица | 
является положительно определенной, то за счет элементов | <” | система 
может потерять устойчивость движения. 

Кососимметричная матрица — диссипации и формирует 
гироскопические силы, а матрица жёсткости | <” | - неконсервативные 
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циркуляционные силы [5-7]. В целом структура сил контактного 
взаимодействия в динамической системе трения такова, что они, за редким 
исключением, приводят к формированию  гироскопических и 
циркуляционных сил. Подчеркнём, что все имеющиеся экспериментальные 
данные по изучению траекторий движения контактируемых поверхностей в 
узлах трения показывают, что при потере устойчивости траектории 
движения индентора всегда остаются круговыми независимо от параметров 
подвески индентора. 

Влияние матриц динамической жёсткости и диссипации процесса 
трения на устойчивость точки равновесия. Проанализируем роль мат- 
риц динамической жёсткости и диссипации в стабилизации и потере устой- 
чивости равновесия динамической системы трения. Отметим главные из- 
вестные особенности связи, формируемой трибосредой [1-3]. 

1. При сближении индентора с контробразцом быстро возрастают 
нормальные составляющие сил, препятствующие сближению поверхностей 
и характеризующие, например, несущую способность подшипника. В неко- 
торых случаях [1-4] имеют место аномальные свойства функции сближе- 
ния, характеризуемые многозначностью силовой функции в зависимости от 
сближения поверхностей. Однако в большинстве случаев эта характеристи- 
ка монотонна и однозначна. 

2. Происходит запаздывание тангенциальных составляющих сил 
при изменении нормальных составляющих. Запаздывание обусловлено 
перестройкой стационарного состояния в трибосреде, которая требует про- 
хождения некоторого пути при движении индентора относительно контроб- 
разца. Имеющиеся экспериментальные данные подтверждают наличие это- 
го запаздывания и его увеличение при уменьшении скорости относительно- 
го скольжения. 

3. При увеличении скорости относительного скольжения при прочих 
неизменных условиях в отдельных диапазонах скоростей происходит 
уменьшение тангенциальной составляющей силы контактного взаимодей- 
СТВИЯ. 


Вышеизложенное позволяет уточнить вектор -— функцию 
ЕСХ,Х,Г) и представить её в следующем виде: 
Е = Е (ХИ); 
(5) 


Е = АЕ (#-1)- Е’ Х). 
С учётом (5) матрицы с; и И; будут иметь вид: 
[а -9А /9ж}; С И й, | 
- (с›-ЮЕ/9ж); с,> — (+ М9Е/9х); [№5+980)/9%,] ^ 
В зависимости от С’ и Г’ значения 9Ё` /дх, монотонно меняются в 


Т,пи Т, 2 
диапазоне [с "< 9Е/9дх |< в а причём 9Е/9ж < 0. 








| . 
Коэффициент 92 Е в зависимости от И, может быть как 
Т 
положительным, так и отрицательным. Обозначим -9Ё /дх, = в и 


Вод о: Тогда система (3) может быть представлена как 
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тх+ МОХ сх И Х+ с®х= 0, (6) 
во в; [^. ы Ин ст] (с) [< у ат ; 16 ы фт 
где = = : 
[%.- Ин]; (в,-№ ' [2+ Ик]; с, 


симметричные части матрицы с; и й; С учётом реакции со стороны 
трибосопряжения (они определяют потенциальные свойства системы); 
. 1 (Г) . - № (Т) 
0; Уи Й и 0; 5 КС 


ИО = = 
: (г). ; ‚ т). 

С ИИ Е. 0 р 0 
кососимметричные составляющие матриц с; и И; ‚ определяющие 


т 
гироскопические и циркуляционные силы. Параметр ср определяется 


градиентом функции сближения в точке равновесия, то есть характеризует 
нормальную контактную жёсткость трибосопряжения. Это обобщение 
понятия контактной жёсткости, принятого в [8]. Коэффициент А можно 
интерпретировать как динамический коэффициент трения, 
рассматриваемый в вариациях относительно стационарного состояния. В 


[1-3] показано, что он может принимать, в зависимости от сближения 
поверхностей, существенно различные значения, даже близкие к единице. 


Запаздывающий аргумент Т также зависит от Х’*. Принципиально, что 

медленно изменяющейся траектории Х*’ соответствует траектория 
и * * * 

параметров с\1’(Х”), АХ") ит(Х°). 


Анализ матриц Й“°”, И, с< и с позволяет определить 


некоторые механизмы потери устойчивости точки равновесия системы и её 
структурных преобразований. 
Первый случай. Рассмотрим влияние асимметрии позиционных сил 


в случае >20 и #->0. Тода #9530, АО=ИЕЛ, 
и вместо (6) имеем: 

тх+ Их+ сОх+ с®х= 0. (7) 

Мы видим, что позиционная связь, формируемая процессом трения, 

вызывает изменение потенциальных свойств динамической системы 

(определяются матрицей с“) и приводит к образованию циркуляционных 


сил (определяются матрицей =: >), которые могут влиять на устойчивость 


точки равновесия. Изменение положительно определённой матрицы Й 
может лишь преобразовать устойчивую по Ляпунову систему в 
асимптотически устойчивую, если матрица жёсткости симметрична и 
положительно определена. Стабилизировать систему, неустойчивую по 
позиционным силам, с помощью изменения матрицы А невозможно. 

Из (7) получаем характеристический полином системы: 


А, (р)=А(р)+ с тр” + (> - АЙ, р+ [=> ы Ас, , (8) 


где А(р) - характеристический полином асимптотически устойчивой 
системы без трения (левая часть (2)). 
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Анализ (8) показывает, что если исходная динамическая структура 
подвески индентора является ортогональной', то циркуляционные силы не 
могут привести к потере устойчивости. Действительно, при 


ср = И, = 0 полином (8) преобразуется В полином 
А. [р] = [ тр? + Ире с [ тр? + №,р+с,.|, корни которого — имеют 


отрицательные вещественные части. Очевидно, что это утверждение 
можно распространить и на случай, когда механическая часть системы без 
трения представлена в виде № -мерной динамической структуры, 
обладающей ортогональными динамическими свойствами. Таким образом, 
влияние циркуляционных сил на устойчивость точки равновесия возможно 
лишь в тех случаях, когда связь между координатами формирует 
замкнутые динамические структуры, условия самовозбуждения в которых 


меняются при изменении параметров матрицы 5” : 


Вначале положим, что К = соп5Ё. Тогда можно вычислить 


предельные значения в и условия, при которых равновесие системы 


является асимптотически устойчивым. Для этого удобно воспользоваться 
частотным критерием устойчивости Михайлова (рис. 2). В зависимости от 
параметров системы, прежде всего недиагональных элементов матриц 


жёсткости и диссипации подвески индентора, а также динамического 


коэффициента трения, при варьировании с“\” возможна потеря 


устойчивости равновесия системы. 





Рис.2. Схема преобразования годографа Михайлова за счёт потенциальных 
и циркуляционных сил, формируемых в трибосопряжении 


Мы видим, что неконсервативные (циркуляционные) силы могут 
приводить к потере устойчивости точки равновесия системы. 
Второй случай. Рассмотрим влияние на устойчивость равновесия 


запаздывающего аргумента 7, изменяющего матрицу диссипации р; и 
формирующего гироскопические силы за счёт кососимметрической части 





' Здесь и далее под ортогональной динамической структурой будем понимать та- 
кую, для которой все недиагональные элемента матриц С и Й равны нулю. 
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р. Как и в первом случае, если И является диагональной, то 
запаздывающий аргумент не влияет на устойчивость точки равновесия 
системы. В общем же случае структура матриц й и И“ в (6) 
показывает, что запаздывающий аргумент Т усиливает положительную 
определённость матрицы АГ“ и дополнительно приводит к формированию 


гироскопических сил, также увеличивающих запас устойчивости системы. 
Для доказательства этого утверждения рассмотрим характеристический 


полином системы (6) в предположении, что Жо = 0. Тогда с учётом (8) 
имеем: 

А : (р)=А,(р)+т Кс р й, р + С.1], (.9) 
где А; (р) - то же, что ив (8). 

Пусть полином А;(р) соответствует системе, устойчивой по 
Ляпунову. Тогда годограф Л, (1@) проходит через начало координат 
(рис.З,а). Каждая точка годографа ЛД’ (/@0) смещается влево на 
_ ФАС) й. ? и вверх на величину ФАСО) с> 107 ‚ То есть во внешнюю 
сторону по отношению к годографу А,;(р), соответствующему 
устойчивой по Ляпунову системе. Следовательно, годограф й (719) 
будет отвечать асимптотически устойчивой системе. Более того, если 
годограф Л’, (70) соответствует неустойчивой системе, то по мере 


увеличения Т система может приобрести асимптотическую устойчивость 
(рис.3,6). 





Рис.3. Влияние на годограф Михайлова запаздывающего аргумента 7 


Третий случай. Рассмотрим влияние на устойчивость параметра #5”, 


который связан с кинетической характеристикой процесса трения [4]. К 
аналогичным эффектам приводит и запаздывание нормальных составляю- 
щих сил по отношению к нормальным смещениям индентора. Анализ мат- 


й й0> 
риц Я; ис; в (6) показывает, что параметр /1>> влияет на симметрич- 
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(с) й92 
ную составляющую й . Более того, при увеличении 22 ВОЗМОЖНО 
У ‚ 


преобразование матрицы и" в отрицательно определённую. В этом слу- 


чае в системе трения могут формироваться ускоряющие силы, при всех 
условиях вызывающие потерю устойчивости точки равновесия системы. 
Выводы. 

1. При взаимодействии механической системы с трибосредой 
объединенная динамическая система обладает принципиально иными 
динамическими свойствами по сравнению с исходной системой без трения. 
При этом в контакте в зоне трения формируются диссипативные, 
ускоряющие, гироскопические, потенциальные и циркуляционные силы. 
Характерно, что диссипативные и потенциальные силы зависят как от 
параметров связи, формируемой узлом трения, так и от параметров 
подвески индентора. Что касается гироскопических и циркуляционных сил, 
то они не зависят от параметров подвески индентора. 

Механизм изменения потенциальных и формирования циркуляцион- 
ных сил связан с отсутствием изменения силовой реакции со стороны три- 
босопряжения на тангенциальное смещение контактируемых поверхностей. 
При определённом соотношении нормальных и тангенциальных составляю- 
щих сил контактного взаимодействия динамическая система за счёт фор- 
мирования циркуляционных сил может потерять устойчивость. Возмож- 
ность потери устойчивости зависит также от свойств подвески индентора. 
Если матрицы жёсткости и диссипации подвески индентора являются диа- 
гональными, то циркуляционные силы не влияют на устойчивость равнове- 
сия. 

Механизм изменения диссипативных и формирования гироскопиче- 
ских сил обусловлен влиянием запаздывания изменения тангенциальных 
составляющих сил по отношению к вариациям нормальных составляющих. 
Это запаздывание всегда направлено на стабилизацию точки равновесия 
системы. 

Механизм формирования ускоряющих сил обусловлен кинетической 
характеристикой процесса трения. Ускоряющие силы могут формироваться 
и за счёт запаздывания нормальных к контактируемой поверхности сил 
при вариациях смещения поверхности в этом же направлении. 

2. Структура формируемых сил такова, что при движении поверхно- 
стей в вариациях относительно точки равновесия по замкнутому циклу со- 
вершается работа не только за счёт сил диссипации, но и за счёт сил упру- 
гости. В связи с этим структура работы, совершаемой силами, обусловлен- 
ными колебаниями в вариациях относительно стационарной траектории, 
достаточно сложная, и это обстоятельство необходимо учитывать при изуче- 
нии эволюционных преобразований в динамической системе трения. 

Выполненный анализ относится к малым колебаниям относительно 
точки равновесия, поэтому при потере устойчивости в системе необходимо 
анализировать нелинейные связи, формируемые процессом в вариациях 
относительно точки равновесия, и образующиеся при этом многообразия в 
пространстве состояния системы. 
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